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MEAN SEA LEVEL DETERMINATION FROM 
SATELLITE ALTIMETRY 


W. D. Kahn 

Goddard Space Flight Center 

B. B. Agrawal 
R. D. Brown 

Computer Sciences Corporation 
ABSTRACT 

The primary experiment on the Geodynamics Experimental Ocean 
Satellite-3 (GEOS-3) is the radar altimeter. This experiment's 
major objective is to demonstrate the utility of measuring the 
geometry of the ocean surface, i.e. the geoid. Results obtained 
from this experiment so far indicate that the planned objectives 
of measuring the topography of the ocean surface with an absolute 
accuracy of ±5 meters can be met and perhaps exceeded. 

The GEOS-3 satellite altimeter measurements have an instru- 
ment precision in the range of ±25 cm to ±50 cm when the 
altimeter is operating in the "short pulse" mode. 
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After one year's operations of the altimeter, data from over 5,000 
altimeter passes have been collected. With the mathematical models 
developed and the altimeter data presently available, mapping of 
local areas of ocean topography has been realized to the planned 
accuracy levels and better. This paper presents the basic data 
processing methods employed and seme interesting results achieved 
with the early data. Plots of mean surface heights as inferred by 
the altimeter measurements are compared with a detailed 1° x 1° 
gravimetric geoid. 
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MEAN SEA LEVEL DETERMINATION FROM 
SATELLITE ALTIMETRY 


INTRODUCTION 

The GEOS-3 altimeter experiment objectives are (1) to determine the feasibility 
and utility of a spacebome altimeter to map the ocean geoid with absolute height- 
resolution of ±5 meters and a relative height precision of ±1 to ±2 meters, 

(2) to determine the feasibility of measuring the deflections of the vertical at 
sea, (3) to determine the feasibility of measuring wave height and (4) to contribute 
to technology leading to a future operational satellite altimeter system with a 

10 cm measurement capability. This paper addresses objective 1 above. 

Satellite tracking contributions to the determination of the current ocean geoid 
have a spatial resolution corresponding to a half wavelength of approximately 

11 to 15 degrees (1250 to 1650 km) . Although the addition of surface gravity per- 
mits local representations with finer resolution, in the l-to-5 degree (110-to- 
550 km) range these data are available over only about 50 percent of the ocean 
surface. 

The GEOS-3 altimeter has demonstrated the capability to determine the fine 
structure of the mean sea surface. When this instrument is operating in the 
short pulse mode a measurement is produced every 4 km along the satellite 
subtrack. With this resolution capability it will be possible to describe sea 



surface topography to a detail of less than 1°, depending on the degree of data 
smoothing used and the spacing between subtracks. 

The GEOS-3 altimeter system therefore fills the gap in data coverage and also 
provides valuable independent knowledge where data now exist. 

When the altimeter data are calibrated and corrected (e. g. , for sea state, 
ocean tides, and other effects) the data constitute measures of the distance 
between the GEOS-3 satellite and the mean sea level, (MSL) or geoid surface. 
Knowledge of the satellite orbital altitude relative to the reference ellipsoid 
will then permit the determination of the geoid. The chief problem encountered 
thus far is the accurate determination of the orbital altitude for GEOS-3. This 
paper will primarily discuss the altimeter measurement geometry and the 
mathematical procedure used to correct orbit errors in the GEOS-3 altimeter 
data. Intercomparisons between processed altimeter sea height profiles and 
independent geoid estimates will be shown. 
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MATHEMATICAL MODELS: 

GEOID PROPERTIES 

The geoid is an imaginary surface, a particular equipotentlal surface of the 
geopotential field. This surface is nearly spherical, with a slight flattening a~ 
the poles and a pear-shaped bulge in the southern hemisphere. Smaller scale 
deviations from the spherical surface can be deduced from the magnitudes of the 
higher order terms in a spherical harmonic model of the geopotential; this 
reveals a nonuniform distribution of surface roughness with surface relief on 
the order of one part in ten thousand. The geoid is approximated closely by a 
real surface, the sea surface, in ocean areas. The sea surface would conform 
exactly to the geoid in the absence of currents, tides, and weather-related 
phenomena (Reference 1). The exact shape of this mean-sea- level (MSL) or 
marine geoid surface is dictated by the spatial variations in the gravitational 
potential of the rotating earth. Its shape is important for this study because 
ideally it is the sea surface that is defined by the altimeter measurements. 

The properties of the MSL surface are represented by the geoid undulation, which 
is the deviation of the geoid from the best fitting ellipsoid as measured along 
the ellipsoid normal. 
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Undulation Amplitude 

Because the magnitude of the geoid surface relief is very small relative to the 
earth's radius, a reference surface that closely approximates the geoid is 
customarily utilized to avoid working with small differences in large numbers. 
This reference surface is usually chosen as the rotationally symmetrical 
ellipsoid that best fits the geoid, in a least-squares sense, to the accuracy 
that the geoid is presently known. This surface plays another role in that it 
represents an equipotential surface of the "normal" potential of the geopotential 
field. The additional factor (the "disturbing" potential) needed to fully de- 
scribe the geopotential is thus made small. 

The geopotential may be described in geocentric earth-fixed coordinates as 

W(X, Y, Z) = U(X, Y, Z) + T(X, Y, Z) (1) 


where 

W » the geopotential 
U = the normal potential 

T = the disturbing potential at the point (X, Y, Z). 

To evaluate the magnitude of geoid undulations, one compares the geoid 
(described by a particular equipotential surface) 
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W(X, Y, Z) = W 0 

with the reference ellipsoid 


of the same potential 


U(X, Y, Z) = U 0 


Uq-Wq 


( 2 ) 

( 3 ) 

I*) 


Following Heiskanen and Moritz (Reference 2) and referring to Figure 1, a 
point P on the geoid is projected along the ellipsoid normal onto point Q on the 
ellipsoid. This distance PQ, denoted by N, is called the geoidal undulation 
at P. 


One may approximate the value of the normal potential U at P by a linear 
relation 

/ du\ a 

Up=Uq+ (i^)q n ' Uq ' TqN (5) 

where Yq is the acceleration of gravity calculated at the ellipoid surface. 

Since by definition 

W p = Up + Tp ( 6 > 

and since 


Wp= W o = U Q 


(7) 



we find by substituting Equations (7) and (6) into (5) that the geoidal undulation 
is given by 
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MEASUREMENT GEOMETRY AND MATHEMATICAL MODELS 


( 8 ) 


The altimeter measurement is nominally the shortest distance between the 
satellite and the sea surface. It is along the normal to the sea surface that 
passes through the satellite. The satellite altimeter antenna is controlled 
to point along the local vertical at satellite altitude. Because equipotential 
surfaces at satellite altitude are even closer approximations to an ellipsoid 
then is the geoid, this local vertical can also be assumed to be the normal to 
an ellipsoidal surface. However, this normal is not generally the same as the 
reference ellipsoid normal due to curvature of the gradient of the normal 
potential. The use of an antenna beamwidth of 1 or 2° insures that the first 
return from the sea surface is within the beamwidth of the antenna. Assuming 
that the measurement is made along the ellipsoid normal at S’, the factors that 
influence the error is the measured altitude are represented in Figure 2. 


Figure 2 shows the beam pattern of a spacecraft-borne altimeter sensor system 
that is oriented in the direction of the gravity acceleration vector (g) at the 
satellite. The height of the satellite is measured in the direction of shortest 
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signal echo, which is very close to the local normal, n, to the geoid (i. e. the 
direction of line segment SP), especially over ocean areas. 

Ideally, the satellite altitude at P should be calculated along the line from P 
to the satellite, of distance h a . However, the location of point P is not known 
Initially, so approximate techniques are used. Although the quality N* is not ex- 
actly equal to N, nor colinear with the SP line segment, the height of the satellite 
above the ellipsoid may be approximated initially as 

h a as h - N' - 6h s - Ah' (9) 

where 

h * magnitude of h, S/C altitude above the reference ellipsoid at S’ 
h a = S/C altitude above the geoid at P 
N’ = geoid height above reference ellipsoid 
Ah $ = deviation of sea surface from geoid 
Ah' = systematic errors in altimeter measurement, 
e. g. , refraction, antenna offset, timing, etc. 

This error sources that affect the altimeter measurement fall into three 
categories: (1) those that result from orbit uncertainly; (2) those that result 
from deviation between the measured surface and the true geoid; and (3) those 
that affect the measurement accuracy itself, i. e. , instrument errors and prop- 
agation effects. The functional dependence of the error sources on the 
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altimeter measurement is given by (Reference 3) 


h a = h(E, t) - N'(G) - 5h s (r, s, 0, X, t) - Ah'(m, 0, X, t) 


( 10 ) 


where 

E = orbit parameters and orbit dependent terms (radiation pressure, 
drag, etc. ) 

t = time of altimeter measurement 
G = vector of geopotential coefficients 
T = vector of tidal error model coefficients 

s’ = vector of local sea-surface biases (i. e. , currents, storm surges, 
wind waves . . . ) 

m = vector of altimeter measurement error coefficients 
0 = latitude 

X = longitude 


or in general 


h a h a (E, G, r, $, m, 0, X, t) 


( 11 ) 


The dominant error source is that due to uncertainty in the orbit. 

The instantaneous sea surface generally is within a meter or two of the geoid 
(Reference 1), but such deviations cannot be neglected if an altimeter accuracy 
of 1 meter or less is desired. Similarly the measurement errors are neglected, 
being of the order of 1 m. 
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Comparing the altimeter measurement, h' with altitude calculated from the 
satellite's orbit h a , the measurement residual Ah a is formed. 

A h a = h '- h a (12) 

The computation of altitude from an orbit based on independent tracking data 
sets has been noted to produce residuals which exhibit systematic offset or 
trends. This systematic offset, due in large measure to orbit uncertainty, is 
reduced by differential correction of the orbit using altimeter data. The re- 
sultant residuals in satellite altitude exhibit almost complete freedom from 
systematic errors. 

The data processing procedure for accomplishing this removal of orbit 
uncertainty effects is as follows. The high rate (10/sec or 100/sec) binary 
format altimeter measurement data is processed in the GEODYN Program, 
a general purpose orbit determination program in use at the NASA Goddard 
Space Flight Center, where a tidal correction is applied and the data is 
converted to geoid heights. The program has been modified to permit the dif- 
ferential correction of orbital elements from altimeter data. This process 
can be used to converge the altimetric geoid heights on a profile through a 
reference geopotential level surface, in this case, the Goddard Earth Model-7 
(GEM- 7) geoid (Reference 4). The adjusted orbit parameters serve merely as ac- 
commodation coefficients in this process and effect a removal of bias , trend, and 
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long wavelength curvature from the relatively short arcs of altimeter data. 
While all six orbit parameters are adjusted, significant changes (greater than 
the standard deviation of the estimate) usually only occur in the orbit semi 
major axis, eccentricity and the true anomaly. Significant changes in the 
remaining Keplerian parameters may also occur for the longer arcs. 

This process minimizes the square of deviations of the altimetric geoid heights 
from the GEM-7 geoid, as is shown in the following development. The ex- 
acted altitude of the spacecraft may be written (neglecting sea surface effects 
and altimeter system errors ) as 

h a (t) = h(e,t)-N'(G) (13) 

where h is the satellite height above the reference ellipsoid as a function of 
the state vector cf the satellite e, and N* is the geoid undulation as a function 
of geopotential parameters G. The actual altitude measurement (t) is re- 
lated to the expected altitude (with suitable differential corrections to the 
parameters) by 


, dh - 3N' - 

h' (t) — h„(t) = — Ae 5 AG + 5 


3e 


3G 


( 14 ) 
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where & is the vector of unmodelled errors in the altimeter measurement. 
Subsequent formation of normal equations solving for corrections to the 
parameters on the right hand side of equation 14 will yield a least squares best 
fit of \ to h; . 

We may rewrite (14) in terms of geoid undulations by substituting (13) into 
(14) and rearranging 

, - - . 3h -* 3N' , 

N (G) = h(e, t) - h' (t) + — Ae AG' + 5 (15) 

3e 3G' 

Due to the near colinearity and similarity of dimension between h a and N', we 
may consider (15) to be an observation equation for geoid undulations. If we 
desire to obtain a least squares fit between these altimetric geoid undulation 
observations and the GEM-7 geoid undulations, we need only substitute N* 
(GEM-7, t) for the left-hand side of (15) and adjust suitable parameters on the 
right hand side. To evaluate the short wavelength differences between the 
altimetric geoid and the GEM-7 geoid, adjustment of £ alone is indicated, since 
these parameters affect only the relatively long wavelength behavior of the 
altimetric geoid data record. Therefore, we adopt the following equation 

N'(GEM-7, t) - h(e, t) + h^(t) = ^ A e + 6 (16) 

3e 

as our observation equation, form the corresponding normal equation, adjust 
the orbit parameters in an iterative differential correction process, and 
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obtain a converged set of altimetric geoid heights 

N(t) = h(e c , t) - ha<t) (17) 

which are a least squares fit to the GEM-7 geoid heights N(GEM-7, t). The 
adjusted orbit parameters e c serve only to effect the fitting process and are 
of no intrinsic significance. Note that equation (16 ) is exactly equivalent to 
equation (14). We have simply transformed our observations from altitudes 
to geoid heights. 

An example of the results of this process is shown in Figure 3 for pass 
number 254. This process removes all bias, linear drift, and long wavelength 
phenomena from the altimetric geoid heights, including possible long wave- 
length geoid information. Since GEM-7 is complete only through degree and 
order 16, with a few terms as high as degree 30, one expects the converged 
geoid heights spectrum to contain geoid information and altimeter noise of 
wavelength shorter than 2500 km. One must be careful in attempting this 
convergence with passes shorter than 2500 km, because aliasing with the long 
wavelength terms in GEM-7 may cause erroneous removal of information of 
shorter wavelengths. 

Implicit in this process is the assumption that the GEM-7 geoid is an accurate 
representation of the long wavelength effects of the geoid, which may or may 
not be true. In any case, it is a simple matter to convert these converged 
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geoid heights to another more accurate analytical model reference surface 
when the need arises. 

Results 

The GEOS-3 altimeter has consistently demonstrated the ability to detect true 
sea surface topography. Figure 4a represents an altimeter pass going from the 
northeast towards the southwest in the western North Atlantic. The satellite 
passes over the well known Puerto Rico trench. This feature is clearly 
indicated in the sea surface height profile as generated from altimeter data 
(Figure 4b). In this same figure very close agreement is seen between the 
GEOS-3 sea surface height profile and the 1° X 1° detailed gravimetric geoid. 

In figure 4c, agreement between the GEOS-3 sea surface height profile and 
bathymetric measurements is also demonstrated. 

Other geographic areas traversed by GEOS-3 also exemplify the aforementioned 
agreement between the altimeter measured sea surface topography and the 
1° X 1° detailed geoid. Figure 5 shows the close comparison between the 
altimeter measured sea surface height profile and the 1° X 1° detailed geoid 
for the mid- Atlantic ridge region. 

Figure 6 shows the fine structure differences between GEOS-3 altimeter 
measured sea surface height profile and the detailed gravimetric geoid. 

Of particular interest is the Aleutian trench area. There is close agreement 
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between the altimeter measured sea surface height profile and the 1° X 1° 
detailed geoid except in the Aleutian trench area where the difference is as 
high as 10 meters. This difference reflects the lack of surface gravity data 
of sufficient accuracy and densification in that region. It suggests that altim- 
etry can make a significent contribution to more accurate geoid/geopotential 
modeling. 

Figure 7 is an altimeter profile of both sea and land topography which clearly 
shows the variation in height between the ocean area and land. Most spectac- 
ular are the signatures of lakes Erie, Huron, and Superior. The variation 
between sea level and the height measured by the GEOS-3 altimeter of these 
Great Lakes closely correlates with topographic survey measurements. 

The large excursions in "sea surface height" are, as one would expect, due 
largely to the changes in received signal strength as the altimeter traverses 
areas of changing surface reflectivity. 

Finally, Figure 8 shows the GEOS-3 altimeter sea surface height profile over 
the Formosa Strait, in the vicinity of the boundary between the Philippine 
and China tectonic plates. The fine structure of the measured altimeter height 
profile does not agree with the detailed geoid for this area, again due to the 
lack of surface gravity data of sufficient accuracy and densification. It 
should be noted, that the altimeter depicts a geoid depression which is not 



18 


associated with a bathymetric feature such as a trench, but which is recognized 
as a subduction zone between the above mentioned plates. 

SUMMARY 

Data from selected GEOS-3 altimeter passes have demonstrated the ability 
to represent precisely the fine structure of the sea surface. Close agreement 
has been shown between the measured sea height profile, the detailed geoid, 
bathymetry, and the structure of tectonic plates. Future application of altim- 
eter measurements will be to determing a detailed global geoid with an absolute 
vertical accuracy of better than 5 meters and a horizontal spatial resolution of 
better than 500 km. 
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Figure 1 Geoid and Reference Ellipsoid 
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Figure 2 Altimeter Measurement Geometry 
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FIGURE 4(«) G EOS-3 ORBIT SUBTRACK FOR 
ALTIMETER PASS 168 


FIGURE 4(c) G EOS-3 SEA SURFACE HEIGHT PROFILE 
COMPARED WITH BATHYMETRY FOR 
PASS 168 


FIGURE 4(b) SEA SURFACE HEIGHT COMPARISON 
WITH 1°x1° GSFC DETAILED GEOID. 
TIME FROM APRIL 21, 197S 
20H 59 M 1SS 


Figure 4 GEOS-3 Sea Surface Height Profile in the Atlantic Area 
Compared with the GSFC 1° x 1° Detailed Geoid 
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FIGURE 5(b) SEA SURFACE HEIGHT COMPARISON 
WITH 1°*1° GSFC DETAILED GEOID. 
TIME FROM MAY 17, 1975 
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Figure 5 GEOS-3 Sea Surface Height Profile in the Atlantic Area 
Compared with the GSFC 1° x 1° Detailed Geoid 






FIGURE 6(a) GEOS3 ORBIT SUBTRACK 
FOR ALTIMETER PASS 
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FIGURE 6(b) SEA SURFACE HEIGHT COMPARISON WITH 1°x1° 
GSFC DETAILED GEOID. TIME FROM MAY 4, 1975 
13H 54M 23S. 


Figure 6 GEOS-3 Sea Surface Height Profile in the Northern Pacific Area 
Compared with the GSFC 1° x 1° Detailed Geoid 
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FIGURE 7(b) TIME FROM JULY 17, 1975 
17H 7M 25S. 


Figure 7 GEOS-3 Sea Surface Height Profile in the Western Atlantic Area 
and Continental United States 
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FIGURE 8(a) GEOS-3 ORBIT SUBTRACK 
FOR ALTIMETER PASS 
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FIGURE 8(b) SEA SURFACE HEIGHT COMPARISON WITH 
1°x1° GSFC DETAILED GEOID. 
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Figure 8 GEOS-3 Sea Surface Height Profile in the Formosa Strait Area 
Compared with the GSFC 1° x 1° Detailed Geoid 




